Epitaxial growth mechanisms of graphene and effects of substrates 
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The growth process of single layer graphene with and without substrate is investigated using 
ab-initio, finite temperature molecular dynamic calculations within density functional theory. An 
understanding of growth mechanisms in the atomic level is provided by exploring the transient 
stages which occur at the growing edges of graphene. These stages are formation and collapse of 
large carbon rings together with the formation and healing of Stone- Wales like pentagon-heptagon 
defects. The activation barriers for the healing of these growth induced defects on various substrates 
are calculated using the climbing nudge elastic band method and compared with that of the Stone- 
Wales defect. It is found that the healing of pentagon-heptagon defects occurring near the edge in 
the course of growth is much easier than that of Stone- Wales defect. The role of the substrate in 
the epitaxial growth and in the healing of defects are also investigated in detail. 

PACS numbers: 68.55.A-, 81.10.Aj, 81.15.-z 



I. INTRODUCTION 

Unusual chemical and physical properties, such as high 
chemical stability, high carrier mobility and high mechan- 
ical strength have made graphend^ an attractive material 
in fundamental and applied science. Pristine graphene, 
in particular, appears to be an important material to 
be used in future high-technology applications.'^ In this 
respect, the production of epitaxial graphene has been 
the motivation of recent experimental and theoretical 
studies-l^ 

The recent studies on graphene production processes 
can be grouped in two classes: one class includes the me- 
chanically exfoliated graphene sheets where th e grap hene 
flakes are peeled from a bulk graphite substrate.'^SISl How- 
ever, that method has the disadvantage of not being 
able to easily control the size and quality of the fabri- 
cated layer. The other class includes the direct growth of 
graphene flakes on substrates.!^ It was shown that epi- 
taxial graphene is a suitable material for nanoscale elec- 
tronic applications by growing ultrathin graphite layer on 
silicon carbide by thermal decomposition.™ Along with 
one layer epitaxial graphene, multi-layer graphene has 
also been studied. 

There are two main events happening during epitaxial 
graphene formation: Nucleation and growth of graphene 
from the nucleated seed. Nucleation on a substrate is 
mostly favored by defects which are actually step edges 
in the atomistic scale. The nucleation process of car- 
bon monomers and dimers on various transition metal 
substrates were investigated in a recent study to distin- 
guish the substrates depending on monomer- or dimer- 
based carbon nucleation.'^ Mechanisms and factors in- 
fluencing graphene growth were also investigated using 
Ru(OOOl), Pt(lll), Ir(lll) and Pd(lll) substrates-'SHIli 
It was found that the orientation of grown graphene 
sheets and carbon concentration at their edge play crucial 
role in the growth process.!^ Graphene was also grown 
on SiC substrates using high temperature sublimation.'^ 



In a recent study, the stability of graphene on nickel 
surfaces and the healing of defects were investigated using 
Monte Carlo simulations with tight-binding potentials.'^ 
Defected graphene flakes containing pentagons, hep- 
tagons, octagon and nonagons were treated at different 
temperatures and the healing effect of increasing temper- 
ature was observed. Also in this same study, nickel was 
proposed as a substrate for better defect healing at low 
temperature, which is consistent with the results of our 
calculations. 

In this paper we investigate the atomistic mecha- 
nisms taking place during the growth of graphene using 
ab-initio finite temperature molecular dynamics (MD) 
method. To clarify the effects of substrates, we con- 
sidered the hypothetical growth of graphene without a 
substrate, as well as the growth on the layered BN 
substrate. li^l Our calculations reveal two major mecha- 
nisms. The first one is the formation of large carbon 
rings at the edges of the growing structure. With the 
inclusion of new carbon atoms, these large rings further 
expand and eventually collapse into smaller structures 
when some critical ring size is reached. The smaller struc- 
tures formed after the collapse are composed of hexagonal 
and defected regions of graphene. The defected regions, 
which form both before and after the ring collapse, gen- 
erally consist of pentagons and heptagons. The second 
major mechanism of the growth process is the healing of 
these defected regions formed at the edges of the growing 
structure. In this respect, we investigate the healing pro- 
cess of defects composed of neighboring pentagons and 
heptagons which is named as pentagon-heptagon (PH) 
defect throughout the manuscript. PH de fect i s similar 
to the well-known Stone- Wales(SW) defect ,'^^□31 but here 
some carbon atoms are two-fold coordinated with unsatu- 
rated sp^-type bonds. The energy barrier for the healing 
of PH defects are lower compared to that of SW due to 
this deficiency. We present also the analysis for the ener- 
getics of healing of SW and PH defects in free standing 
graphene, as well as graphene grown on BN and Ni(lll) 
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substrates using the climbing image nudged elastic band 
(NEB) method:20i 



II. METHOD 

We have performed atomic structure optimizations and 
ab-initio finite temperature molecular dynamics (MD) 
calculations within density functional theory (DFT) us- 
ing VASP softwarePil The relaxed geometries of all 
structures were calculated by spin-polarized plane-wave 
calculations using projector augmented-wave (PAW) 
potentialJ^ within generalized gradient approximation 
(GGA;P including van der Waals corrections^^ The 
Brilliouin zone of the primitive unitcell of graphene 
was sampled by (17x17x1) k-points in Monkhorst-Pack 
scheme which was scaled according to the size of other 
unitcellsP The energy convergence value between two 
consecutive steps was chosen to be 10^^ eV. In ab-initio 
MD calculations the time step was taken 2.5 fs and the 
atomic velocities were renormalized to the temperature 
set at T=1300 K at every 40 time steps. The temperature 
of MD calculations is in compliance with the temperature 
used in chemical vapor deposition (CVD). 



III. MOLECULAR DYNAMICS SIMULATIONS 
OF GROWTH 

A. Growth without Template 

In order to understand the effects of a template surface 
on graphene growth, we first consider a hypothetical situ- 
ation and investigated the growth mechanism of graphene 
without any template (substrate). To bypass the initial 
nucleation process, a graphene flake was fixed in space 
and additional carbon atoms were sent to it in differ- 
ent scenarios, which mimic the growth. Simulations were 
done by letting the carbon atoms to move freely in a cer- 
tain plane, while not allowing the out of plane motions. 
As shown in Fig. [ija) , we start with nanoribbons having 
armchair and zigzag edges. Initially, zigzag nanoribbon 
is composed of 24 atoms forming three zigzag chains in 
the periodic direction, two of which are kept fixed. Arm- 
chair nanoribbon starts with 20 atoms, 12 of which are 
kept fixed. Fixed atoms are delineated inFig.jlJa). After 
running 1 ps of MD simulation two more carbon atoms 
are introduced in both systems. These atoms are first 
positioned in the same plane to the left of the nanorib- 
bons, where the edges are free to move in 2D. Then they 
are moved towards these edges until the distance between 
the newcoming atom and one of the edge atoms decreases 
to 1.3 A. Then the new MD simulation is started and the 
process is repeated consecutively. 

As seen in Fig. [iJa), during the growth simulation the 
orientation of the honeycomb parts are changed from 
zigzag to armchair (left panel) and vice versa (right 




FIG. 1: Snapshots from ab-initio MD simulations of planar 
graphene growth at T=1300 K without a template substrate. 
An initial flake was placed and in each 1 ps MD calculation 
two carbon atoms were sent from the left hand side to mon- 
itor the growth in the indicated direction. Each snapshot 
includes two periodic supercells. (a) Change of the armchair 
edge to zigzag edge and vice versa is shown, (b) Structures 
obtained when simulation of growth presented in (a) is pro- 
ceeded. Formation of big rings and chains were observed, and 
the resulting structures were far away from being a perfect 
graphene layer. Note that defects formed in part (a) are still 
present in part (b). 



panel). Interestingly, in both cases the transition is me- 
diated by similar structures composed of two heptagons 
with one pentagon in the middle. One can attribute 
the defected growth to the absence of a substrate which 
would act as a stencil if the template had a structure 
similar to graphene. Carrying on the growth simula- 
tion of structures presented in Fig. [ija) results in the 
massively defected network of carbon atoms as shown in 
Fig. [ijb). One can identify the big holes surrounded by 
carbon chains and patterns composed of pentagons and 
heptagons. It was observed that the defects that emerged 
at the beginning of the simulation are still present after 
about 40 ps of simulation. This implies that the process 
of growth induces defects mainly composed of pentagons 
and heptagons, which persist due to the absence of a 
healing mechanism. Several scenarios of growth simula- 
tion without a substrate were tested but none of them 
resulted in a reasonably ordered honeycomb structures. 
Especially, those simulations in which carbon atoms were 
allowed to move in all directions resulted in bulk-like 
structures and atomic chains. Although the observed 
chain structures are interesting, growth of regular hon- 
eycomb structure was not observed and we deduced the 
necessity of a template during the growth process to de- 
fine a plane where the graphene like structure sits and 
where the newcoming carbon atoms are landed. 
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FIG. 2: Snapshots from ab-initio MD simulation of epitaxial growth of graphene on a BN substrate. In the ball and stick 
model B, N and C atoms are represented by green, blue and brown balls while only bonds between carbon atoms having 
distance less than 2 A are shown. Each snapshot includes two periodic supercells in the horizontal direction, (a) General 
trends are presented by including final configurations of MD calculations involving 30, 35, 40 and 42 carbon atoms. Some of 
the critical configurations in the evolution of ring collapse and defect healing mechanisms are highlighted by solid and dashed 
lines respectively, (b) Snapshots from the MD simulation of the structure having 42 carbon atoms taken after 1, 7, 14 and 20 
MD steps. Carbon atom migration causing the growth of rings and defect healing can be traced in dotted and dashed circles, 
respectively, (c) Snapshots from the same MD simulation taken after 40, 83, 222 and 290 MD steps. Three subsequent hexagon 
formations are indicated by solid, dashed and dotted circles. 



B. Dynamics of Graphene Growth on BN 
Substrate 

Hexagonal boron-nitride consists of single layers of 
BN in honeycomb structure, which is almost commen- 
surate to graphene. It has been argued that BN lay- 
ers of any thickness can be grown on graphene layers 
and vise versa.^ Because of this reason we have chosen 
BN as a template on which we investigate the growth of 
graphene. Again, to skip the initial nucleation process, a 
graphene flake was placed on BN substrate. Single car- 
bon atoms were released from random positions on top of 
the graphene flake edges and molecular dynamics simu- 
lation were performed for 400 time steps before the next 
atom was sent. By sending the atoms one by one, events 
happening during the growth process was monitored at 
atomistic scale. The snapshots taken from this calcula- 
tion are shown in Fig. [2j Normally the bottom of each 
panel comprises The bottom part of these structures nor- 
mally comprises fixed graphene and BN substrate, which 
are not shown while growth proceeds upwards. 

In Fig. [2|a) we present the general trends observed 
during growth. Each of the four snapshots in this row 
corresponds to the final structures obtained after the MD 
simulations of 30, 35, 40 and 42 carbon atoms on BN re- 



spectively. As seen in the first column of Fig. [2ja) first 
a single carbon atom indicated by a small arrow makes 
bonds with armchair edges and a pentagon structure is 
formed. This stretches the edges and prepares a medium 
for the formation of a neighboring heptagon. Then ring- 
like structures start to grow at the edges as seen in col- 
umn two of Fig. [2]ja). When the ring structure reaches a 
certain size, it collapses forming hexagonal structures at 
the graphene edges as shown in columns three and four. 

The formation and healing of pentagon-heptagon de- 
fects, which is the second major mechanism affecting the 
growth process is presented in Fig. [2]Jb). The snapshots 
correspond to the 1^*, 7*'*, 14*'' and 20*'' MD steps of 
the simulation with 42 carbon atoms. The healing of the 
PH defect is highlighted by dashed circles. As simula- 
tion proceeds to the 20*'' step, the PH defect is totally 
relaxed into two hexagons. Note that, healing of PH de- 
fect at the edge is similar to the healing of SW defect, 
which involves rotation of the middle bond (which is com- 
mon to two adjacent heptagons) by 90°. Here there is, 
however, a crucial difference in the path of healing com- 
pared to that of SW healing, because one of the carbon 
atoms in the pentagon of the PH defect is bonded with 
two adjacent carbon atoms, whereas in the SW defect 
all carbons are bonded to three others. The absence of 
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one of these bonds (or the presence of sp^-type danghng 
bond) decreases the barrier of PH heahng compared to 
the SW case. This issue is revisited in forthcoming de- 
tailed discussions. 

In Fig.[2jb), the dotted circle in the first column marks 
the inclusion of the newly added carbon atom which is 
added from a random position on top of the graphene 
layer. This newly added carbon atom is initially posi- 
tioned on top of another carbon atom of the graphene 
structure. It then migrates to the bridge site and by re- 
placing the bridge bond it increases the size of the carbon 
ring (column 2 to 4). As the ring expands with the inclu- 
sion of this new carbon atom, it reaches the critical size 
after which it collapses. 

Earlier, the zero-temperature, static calculations based 
on the minimization of the total energy predicted that a 
single carbon atom is preferably adsorbed at the bridge 
site on graphene surface with a binding energy of 2.3 
gy jl?! -^Yhile the energy barrier for the migration of sin- 
gle, isolated carbon atom is only 0.37 eV, it is lowered and 
eventually collapsed when another carbon atom is located 
at close proximity. Both first-principles total energy and 
finite temperature MD calculations have shown that ini- 
tially C2 and eventually C„ carbon chains can form per- 
pendicularly attached to graphene surface through inclu- 
sion of single carbon adatom one at a time. The gain of 
energy in the implementation of a single carbon adatom 
is ~ 5 eV.f^^ However, the situation is dramatically differ- 
ent for a graphene sheet having armchair or zigzag edge: 
Single carbon atoms have shown to be attached favor- 
ably to the edge atoms with much larger binding energy 
(7.08 eV for armchair edge and 8.19 eV for zigzag edge). 
Further implementation of carbon adatoms to the edges 
give rise to PH-hke defect structures. Hence, the results 
of earlier static, zero-temperature calculation^^!! confirm 
the findings obtained from the present finite temperature 
MD calculations. 

The process of ring collapse is shown in Fig. I2|c). The 
snapshots presented here correspond to the 40^'', 83'''', 
222*'' and 290*'' MD steps of the simulation with 42 
atoms. Here the ring is composed of 14 carbon atoms 
before the collapse. This is just enough to form three 
hexagons highlighted by solid line in the fourth column 
of Fig. [2]; a). As seen in Fig. [2jc), these three hexagons 
are consecutively formed during the collapse of the ring. 
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FIG. 3: (a) Energetics of the healing path of SW and PH 
defects in graphene for three cases; namely without template, 
graphene on BN and graphene on Ni(lll) substrates. The 
solid red, green and blue lines show the healing path of SW 
defect and associated energy barriers for graphene without 
template, graphene on BN and graphene on Ni(lll) surfaces, 
respectively. Energies of the defected states are set to zero. 
The PH healing barriers are also shown with dashed lines. 
The inset shows the healing of PH defect of graphene on BN 
substrate. In the PH defect one of the carbon atoms forming 
the middle bond is missing the third sp^-like bond, (b) Top 
and side views of SW defect healing on Ni substrate. The 
Ni atoms forming the top, middle and bottom atomic lay- 
ers of the substrate are indicated by numerals 1, 2, and 3, 
respectively. The lateral positions of atoms in these layers 
are indicated by sites 1, 2, and 3. The interaction between 
graphene and Ni(lll) is manifested in the side view of the 
fifth NEB image, where carbon atoms forming the C-C bond 
between two heptagon are pulled down when they are passing 
over site-2 and site-3 of the Ni substrate. 



IV. ENERGETICS OF PH AND SW DEFECT 
HEALINGS 



Having noticed the role of PH defect healing in 
graphene growth, we move on by investigating the ener- 
getics and dynamics of SW and PH defect healing in free 
standing graphene, as well as graphene grown on BN and 
Ni(lll) surfaces. The role and importance of topological 
SW defects on fuUerene growth, isomerization and plas- 
ticity of various carbon nano structures had been previ- 
ously reported.lSSlEfilThe presence of SW defects have also 



been observed along with the unique out of plane wave- 
like defect in graphene and other sp^-bonded materials .'^ll 
SW defects not only change the geometrical properties 
of graphene but also influence its electronic and chemical 
properties, such as band structure, reactivity and carrier 
transport. 

Here, we calculate the energy barrier confronted dur- 
ing the healing process using the climbing image NEB 
method. ^'^ Structures involved in this calculations are 
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composed of armchair graphene nanoribbons with fixed 
edges and defects in the middle. The height between 
graphene nanoribbon edges and substrate underneath is 
set to the optimized value found in the case of infinite 
graphene sheet on infinite substrate. Also the relative 
position is derived in similar way. The optimum config- 
uration of graphene on BN substrate is achieved when 
carbon atoms of one graphene sublattice are placed on 
top of boron atoms of the underlying BN layer. In case 
of the Ni(lll) substrate, graphene structure is oriented 
in such a way that nickel atoms at the top layer of the 
substrate are under the center of the bridge bonds of 
graphene. We first calculate the ground state configu- 
ration of completely defected and healed states. As an 
initial guess of a healing path, we choose a straight line 
connecting these defected and healed states via linear 
interpolation. We choose 11 NEB images where first (de- 
fected) and eleventh (healed) are not changed while other 
nine images are varied until the optimum healing path is 
found. The fifth image was chosen as the climbing im- 
age which converges to the saddle point. Results of these 
calculations are outlined in Fig. [3j 

In Fig.|3][a), the calculated healing paths of SW defect 
in graphene and associated barriers are shown for three 
cases; namely for graphene without substrate, graphene 
on BN and on Ni(lll) surfaces. Here the energy barri- 
ers along the healing paths of SW defects are found to 
be 3.74 eV, 3.57 eV and 2.96 eV for free-standing, BN 
and Ni substrate cases, respectively. The energy barrier 
is significantly lowered by the substrates. The effects of 
substrates are proportional to their interaction energy 
with graphene structure. In this respect, the binding en- 
ergies of graphene on BN and Ni substrates are found 
to be 0.13 eV and 0.41 eV per two C atoms, respec- 
tively. How the substrate can lower the barrier energy is 
explained by the top and side views of atomic configu- 
ration of SW defected graphene on Ni(lll) substrate in 
Fig. [sjb). Three layers of Ni(lll) forming an A,B, and 
C stacking of the fee structure are indicated by numerals 
1, 2 and 3 starting from the top layer in the side view in 
Fig. [3] (b). Lateral positions of the atoms of these lay- 
ers are indicated by sites 1, 2, and 3 in top view in the 
same figure. Site-2 and site-3 are energetically favorable 
sites for graphene atoms above Ni(lll), since the bind- 
ing energies of a single carbon atom on site-2 and site-3 
is more favorable compared to that of site-1 by 2.48 eV 
and 2.46 eV, respectively. Here during the healing pro- 
cess of SW defect the energy barrier is lowered because of 
two reasons. The first reason is that, the C atoms which 
form the defect are pulled by site-2 and site-3 of the Ni 
substrate as shown in side view in Fig.lstb). This pulling 



is in the same direction with healing path. The second 
reason is that, by pulling the carbon atoms out of plane, 
Ni substrate increases the distance between these atoms 
and thereby decreases stress in the carbon-carbon bonds 
during the healing. 

The healing paths and energy barriers of PH defect in 
graphene are also shown by the dotted curves in Fig. |3ja) . 
The inset shows the atomic configuration of PH defected 
graphene on layered BN substrate. Unlike SW defect, 
here PH defect has one carbon atom with a sp^-type dan- 
gling bond. The energy to be gained from the saturation 
of this dangling bond by forming a bond with a nearest C 
atom of the heptagon ring becomes the deriving force for 
the healing. As a result, the energy barrier to heal the 
PH defect is lowered dramatically by ^ 2 eV as shown in 
Fig. ^a). Our argument is justified by the fact that the 
barrier lowering occurs also for the healing of PH defect 
in free standing graphene. However, contrary to the ef- 
fect of substrate in the healing of SW defect above, the 
barrier lowering effect of Ni(lll) substrate is weaker than 
that of BN substrate. This is due to the interaction be- 
tween the substrate and sp^-type dangling bond; namely 
the stronger the interaction with substrate, the lower is 
the gain of energy upon saturation of the sp^-type dan- 
gling bond of a two-fold coordinated carbon atom. As a 
result, highest barrier lowering takes place in the healing 
of PH defect in free standing graphene. 



V. CONCLUSION 

Graphene growth and energy barrier calculations of de- 
fect healing were investigated using ab-initio MD calcula- 
tions. It was found that there are two mechanisms which 
play crucial roles in the growth of graphene. First mech- 
anism is the formation of large carbon rings at the edges 
which eventually collapse to form honeycomb structure 
with defects. This collapse is found to be initiated by 
the new coming carbon atoms which replace one of the 
bonds in the ring, and expands it until the critical size 
is reached. Second mechanism is the formation of PH 
defects near the edge and their healing. We have shown 
that the energy barrier needed to overcome during heal- 
ing of the PH defects are much lower than that of the SW 
defects. We have shown that the presence of a BN or Ni 
substrate have crucial effect on growth. These substrates 
guide the formation of honeycomb structures from car- 
bon rings and enable the healing of specific defects as 
growth proceeds. 
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